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We present new spectroscopic and photometric observations of the HAT-P-1 
planetary system. Spectra obtained during three transits exhibit the Rossiter- 
McLaughlin effect, allowing us to measure the angle between the sky projections 
of the stellar spin axis and orbit normal, A = 3?7 ± 2?1. The small value of A for 
this and other systems suggests that the dominant planet migration mechanism 
preserves spin-orbit alignment. Using two new transit light curves, we refine the 
transit ephemeris and reduce the uncertainty in the orbital period by an order 
of magnitude. We find a upper limit on the orbital eccentricity of 0.067, with 
99% confidence, by combining our new radial-velocity measurements with those 
obtained previously. 

Subject headings: techniques: radial velocities — planetary systems: formation — 
stars: individual (HAT-P-1, ADS 16402A) 



Introduction 



Prior to 1995, it was expected that Jovian planets around other stars would inhabit 
wide, circular orbits similar to the Solar System gas giants. It was therefore a surprise when 
the first exoplan et was discovered with a minimum mass of 0.468 Mj up and a semimajor axis 
of only 0.05 AU (IMayor fc Queloall995l ). Since then, 85 "hot Jupiters" — Jovian planets with 



perio ds <10 days — have been detected around Sun-like stars (IButler et al.ll2006l ; iTorres et al. 



20081 ). It is unlikely that these planets formed in situ due t o the low surfac e densities and 
high temperatures of the inner regions of circumstellar disks ( Lin et allll996l ). A more likely 
scenario is that these massive planets formed at a distance of several astronomical units, and 
then migrated inward to their current locations. 

Theories for the inward migration of planets can be divided into two broad categories. 
The first category involves tidal interactions between the planet and a remaining gaseous disk 
(ILin et al.l Il996l : iMoorhead fc Adamsl 120081 ) . The sec ond category involves few-body grav- 



itatio nal dynamics, such as p lanet-planet scattering 



20081 ). dynamical relaxation (jPapaloizou &: Terquem 



Kozai cycles accompanied by tidal fric tion (IHolman et al. 



2001 



Rasio fc Fordlll996l: IChatteriee et al. 



Adams fc Laughlin 



1997 



200 



^) , and 



Fabrycky fc Tremaind 12007 



Wu et al.l 120071 : iNagasawa et al.ll2008l ). One possible way to distinguish between these cat- 
egories is to examine the present-day alignment between the stellar rotation axis and the 
planetary orbital axis. Assuming that these axes were initially well al igned, disk-planet tidal 
interactions would preserve this close alignment (IWard &: Hahnlll994j ). while the second cat- 
egory of theories w o uld a t least occasionally result in large misalignments. For example, 
Adams &: Laughlin! (120031 ) predict a final inclination distribution for dynamically relaxed 
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planet ary s ystems that pea ks near 20° and extends to 85°. Likewise, iFabrycky fc Tremaine 
( 120071 ) and IWu et al.l (120071 ) simulated systems of planets with randomly aligned outer com- 
panions and found that the Kozai interaction resulted in a wide distribution of final orbital 
inclinations for the in ner planet, with r etrogr ade orbits (A > 90°) not uncommon. Similar 
results were found by iNagasawa et al.l (120081 ). for the case in which Kozai oscillations are 
caused by an outer planet, rather than a companion star. 

Spin-orbit alignment can be measured by taking advantage of the Rossiter-McLaughlin 
(RM) effect that occurs during a planetary transit. As the planet blocks portions of the 
rotating stellar surface, the star's rotational broadening kernel becomes asymmetric and its 
spectrum appears to be anomalously Doppler-shifted. The RM effect has previously been 



observed and modeled for eight transiting planetary syst e ms jQueloz et al. 2000 ' L Winn et al. 



2005 



2008 



200 



6. 



2007c 



Winn et al. 



20081 ). In this work, we add HAT-P-1 to this sample. 



Wolf et al.ll200~7l ; iNarita et al.ll2007al lbl: iBouchv et al.ll2008l : lLoeillet et al. 



HAT-P-1 (ADS 16402B) is a member of a GOV/GOV visual binar y and harbor s a sho rt- 
period, Jovian planet. The transits of HAT-P-lb were discoverd by iBakos et al.l (120071 ) as 
part of the Hungarian-made Automated Telescope Network (HATNet). The planet has a 
4.465 day orbital p eriod, a mass of 0.53 Mj up , a radius Rp = 1.20 i?j up (IBakos et al.l 120071 : 
Winn et al.l l2007al ). We have monitored HAT-P-1 using precise radial velocity (RV) and 
photometric measurements made both in and out of transit in order to measure the RM 
effect and improve the precision with which the system's orbital parameters are known. In 
the following section we describe our observations and data reduction procedures. In §3] we 
present the transit model that we fit to our observations, and in § H] we present our results, 
and we conclude in § [5] with a brief discussion. 



Observations and Data Reduction 



2.1. Radial Velocity Measurements 



We observed t he optical spectr um of HAT-P-1 using the High Resolution Echelle Spec- 
trometer (HIRES, IVogt et al.l 11994) on the Keck I 10m telescope and the High Dispersion 



Spectrograph (HDS. lNoguchi et al. 



20021 ) on the Subaru 8m telescope. We set up the HIRES 



spectrometer i n the same manner that has been used consistently for the California-Carnegie 
planet search (IButler et al.l 1 19961 ; iMarcy et al.l 120051 ) . This is also t h e sam e setup that was 
used to gather the 9 Keck/HIRES spectra reported by IBakos et al.l (120071 ) . Specifically, we 
employed the red cross-disperser and used the I 2 absorption cell to calibrate the instrumental 
response and the wavelength scale. The slit width was set by the Of! 85 B5 decker, and the 
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typical exposure times ranged from 3-5 min, giving a resolution of about 60,000 at 5500A 
and a signal-to-noise ratio (SNR) of approximately 120 pixel -1 . We gathered 3 spectra on 
several nights when transits were not occurring, in order to refine the parameters of the 
spectroscopic orbit. In addition we gathered a dense time series of spectra on each of two 
nights, UT 2007 July 6 and UT 2007 September 2, when transits were predicted to occur. On 
each night we attempted to observe the star for many hours bracketing the predicted transit 
midpoint, but there were interruptions due to clouds and pointing failures. However, both 
nights of data provide good phase coverage of the entire transit event. In total we obtained 
79 new Keck/HIRES spectra, of which 49 were observed while a transit was happening. 

For our Subaru/HDS spectra we employed the standard I2a setup of the HDS, covering 
the wavelength range 4940-6180A with the I 2 absorption cell. The slit width of 0"8 yielded 
a spectral resolution of ~45,000. The typical exposure time was 10 min resulting in a SNR 
of 120 pixel -1 . Our Subaru observations took place on 3 different nights spread out over 2 
months. Two of the nights were not transit nights; we gathered 8 spectra on those nights in 
order to refine the parameters of the spectroscopic orbit. The last night, UT 2007 Septem- 
ber 20, was a transit night, and we gathered 25 spectra over 7.3 hr bracketing the predicted 
transit midpoint, of which 16 were gathered during the transit. 



We performed the Doppler analysis with the algorithm of iButler et al.l (|1996l). For th e 
Subaru data we used a version of this algorithm customized for HDS by ISato et al.l (120021 ). 
We estimated the measurement error in the Doppler shift derived from a given spectrum 
based on the weighted standard deviation of the mean among the solutions for individual 
2 A spectral segments. The typical measurement error was 3 m s" 1 for the Keck data and 
7 m s -1 for the Subaru data. The data are given in Table [1] and plotted in Figs. [1] and HI 
Also given in that table, and shown in those figures, are dat a based on the 9 K eck/HIRES 
spectra and 4 Subaru/HDS spectra obtained previously by iBakos et al.l (120071 ). We note 
that the RV timestamps reported by Bakos et al. (2007) are incorrect. They were said to 
be Heliocentric Julian dates, but they are actually Julian dates. We provide the corrected 
dates in Table [H 



2.2. Photometric Measurements 

We obtained photometric measurements of HAT-P-1 during the transit of UT 2007 Oct 8 
using the Nickel lm telescope at Lick Observatory on Mount Hamilton, California. We 
used the Nickel Direct Imaging Camera, which is a thinned Loral 2048 2 CCD with a 6.3' 
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square field of view0. We observed through a Gunn Z filter, and used 2x2 binning for an 
effective pixel scale of 0'/37 pixel -1 . The exposure times varied depending upon conditions 
but were typically 10-12 s, with a readout and setup time between exposures of 34 s. The 
conditions were clear for most of the transit with ~ l'.'O seeing. However, observations 
during ingress were partially obscured by clouds and the data from that time period proved 
to be significantly noisier than the rest; we have excluded those data from our analysis. We 
determined the instrumental magnitude of HAT-P-1 relative to two comparison stars using 
an aperture with an 11 pixel radius and a sky background annulus extending from 15 to 18 
pixels. 

We observed the transit of UT 2007 September 20 with the MAGNUM 2m telescope on 
Haleakala, in Hawaii fcobavashi et al.lll998l : lYoshiill2002l : IVoshii et alJ2003h . The MAGNUM 
photometric observations were conducted on the same night as the Subaru/HDS transit 
observations described in § 12.11 We employed the Multicolor Imaging Photometer (MIP), 
using a 1024 2 SITe CCD with a pixel scale of / /277 pixel -1 . The camera's field of view is 
1 '5, which is much smaller than the field of view of the detector. During each exposure, the 
field was shifted on the detector along a 3 x 3 grid, which allowed us to increase the duty 
cycle since the chip was read out only once for every 9 exposures. Observations were made 
through a Johnson V-band filter, and the exposure times were 10 s, with 40 s per exposure 
for readout and se t up. T he MIP images were reduced with the standard pipeline described 
by lMinezaki et al.l (120041 1 . We determined the instrumental magnitude of HAT-P-1 relative 
to its visual binary companion, ADS 16402A, using an aperture radius of 15 pixels, and 
estimated the sky background level with an annulus from 20 to 25 pixels. 

The photometric data are given in Table [2] and plotted in Fig. [2j In the final light 
curves, the root-mean-squared (rms) relative flux, outside of transits, is 0.0019 for the Nickel 
data and 0.0016 for the MAGNUM data. 



3. The Model 



3.1. An Updated Ephemeris 



The extended time baseline of our new photometric measurements allows us to refine 
the transit ephemeris. We first computed midtransit times from the light curves using the 
method described by I Winn et al.l (j2007al ). In particular, to assign proper weights to the 



1 This is the same camera used by IWinn et al.l (|2007aT ) . which they mistakenly described as a 2048 2 
Lawrence Labs CCD with a 6.'1 x 6.'1 field of view. 
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photometric data during the light curve fitting procedure, we applied a correction to the 
uncertainties to take into account time-correlated noise ("red noi se" ) , which was dete rmined 
by examining the rms residuals in time-averaged light curves (see lWinn et al.l (l2007al )). This 
resulted in a factor-of-two increase in the error bars, relative to a situation in which correlated 
noise is ignored. 



We employed the sa me modeling procedure described in detail by iHolman et al.l (120061 ) 
and I Winn et al.l (j2007bl ) and summarized as follows. We modeled the path of the planet 
across the stellar disk using a parameterized model based on a planet and star in a Keplerian 
orbit abou t their center of mass. We fitted the photometric observations using the analytic 
formulas of lMandel &: Agoll (120021 ) and a quadratic lim b-darkening: law with fixed coefficients 
a based on the tabulated calculations of Claret ( 20041 )1^1 The free parameters were the scaled 



stellar radius R*/a (where a is the semimajor axis), the planet-to-star radius ratio R p /R+, 
the orbital inclination i; and for each light curve, the midtransit time T c , the mean out-of- 
transit flux, and a time gradient of the out-of-transit flux (to account for some systematic 
errors in the photometry). The model fit was carried out using a Markov Chain Monte 
Carlo (MCMC) algorithm with 10 6 links, in which a single randomly-chosen parameter was 
perturbed at each link, with a perturbation size tuned such that ~40% of the jumps were 
executed. The mean values and standard deviations of the posterior probability distributions 
(which were nearly Gaussian in this case) were adopted as the "best-fit" parameters and 
uncertainties. 

We fit a linear ephemeris to al l of the time s listed in Table [31 which includes the new 
transit ti mes and those measu red by I Winn et al.l (j2007al ). We found that two of the entries in 
Table 3 of I Winn et al.l (j2007al ) were incorrect: the first time was wrong because the data had 
not been normalized correctly, and the sixth time was too small by one period because of a 
rounding error in the computer code that generated the table. The corrected times are given 
in TableEJ A linear fit to the transit times had x 2 — 10.2 and 8 degrees of freedom, indicating 
an acceptable fit. The fit residuals are plotted in Figure [3] and the updated ephemeris is 
given in Table HI The uncertainty in our updated period is about 10 times smaller than 
the previous estimate. In our subsequent analysis we fix the period at this value, as the 
uncertainty is negligible for our purposes. 



2 For the Z band, the coefficients were az = 0.18 and bz = 0.34. For the V band, the coefficients were 
ay = 0.40 and b v = 0.32. 
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3.2. The Orbital Eccentricity 



(IBakos et al.l 120071 ) reported a tentative detection of a nonzero orbital eccentricity, e = 
0.09 ± 0.02, based on an analysis of 13 RV measurements. With our expanded RV data 
set, we can check on this tentative detection. We modeled our radial velocity measurements 
using a Keplerian orbit with 6 free parameters: the velocity semiamplitude K, the orbital 
eccentricity e, the argument of pericenter a;, and an additive velocity for each of the 3 velocity 
groups (our Keck v elocties and those o f IBakos et al.l (120071 ); our Subaru velocities; and the 
Subaru velocities of IBakos et al.l (120071 )). The time of transit, T c , and the orbital period P 
were held fixed at the values determined from the photometric data. 

To avoid complications at this stage due to the RM effect, we fitted only those 43 veloc- 
ities that were gathered well outside of transits. Specifically we excluded all velocities that 
were measured within a 4 hr window centered on the calculated midtransit time (the actual 
transit duration is 2.8 hr). To assign proper weights to the RV measurements we needed 



(Saar et al. 


1998; 


Wrieht 


2005) 



We found it necessary to add (in quadrature) 3.7 m s 1 to the measurement errors in order 
to obta i n a x f, of unity. T his jitter estimate i s consistent with the 3.4 m s _1 predicted by 



Wrightl (120051 ) and used by lBakos et al.l (120071 ). In the modeling procedure described in the 



rest of this section, we used the augmented error bars, while Table [T] gives only the inter- 
nal measurement uncerta inties. Unlike our previous analyses of HD 189733 and HD 147506 



( IWinn et al.ll2006l . l2007d ). we found no evidence for a higher night-to-night jitter compared 
to the intra-night jitter. We therefore did not modify the error bars any further than the 
quadrature addition of our jitter estimate. 

We employed an MCMC fitting algor ithm using 10 6 st eps and perturbation sizes re- 
sulting in a 30-50% acceptance rate (e.g. IWinn et al.l 120051 ) . The orthogonal parameters 
describing the eccentricity e and argument of periastron uo were ecoscu = 0.003 ± 0.013, 
esincu = 0.004 ± 0.025. The orbital eccentricity of the HAT-P-1 system was found to be 
smaller than 0.067 with 99% confidence. This is consistent with the theoretical expectation 
that the orbit should have circularized due to tidal friction. The circu larization timescale 
is ~0.23 Gyr assumin g a tidal quality fa ctor of 10 6 (IBakos et al.l 120071 ). and the estimated 
stellar age is 2.7 Gyr (ITorres et al.l 120081 ) . In what follows we assume e = exactly. 
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3.3. Joint Analysis of Radial Velocities and Photometry 

To determine the projected spin-orbit angle and its uncertainty, we simultaneously fitted 
a parametric model to the RV data as well as a composite transit light curve, ge nerated from 



all of t he Z and z photometric data at our disposal, from this work and from IWinn et al. 



( I2007al ). The composite light curve has 1 minute bins, and an out-of-transit rms of 0.00057. 
It is shown in Fig. H] along with the transit RVs. Although our main interest is in the 
spin-orbit parameters, which are largely determined by the transit RV data, we included 
the photometric data in the fit as a convenient way to account for the uncertainties in the 
photometric parameters and their covariances with the spin-orbit parameters, although in 
practice these covariances proved to be small. 

The aspects of the model that attempt to fit the photometry, and the orbital Doppler 
data, have already been described. To calculate the radial velocity during transits, we must 
calibrate the relationship between the "anomalous Doppler shift" that is returned by our 
code for measuring Doppler shifts, and the physical p arameters and con figuration of the star 



and planet. For this purpose we used the technique of IWinn et al.l (120051 ). in which simulated 
stellar spectra are created that exhibit the RM effect, and then these spectra are analyzed 
with the same Doppler-measuring code that is used on actual data. Such simulations are 
needed because the algorithm for measuring Doppler shifts involves fitting for parameters 
that are intended to describe the time- variable instrumental profile of the spectrograph, and 
these parameters may interact with the spectral disortion of the RM effect in ways that are 
hard to predict. 

In our simulations the physical configuration of the planet and star is characterized by 
the transit flux decrement, e, and the velocity of the occulted portion of the stellar disk 
(the "sub-planet velocity"), denoted by v p . We created simulated spectra with the same 
data format and noise characteristics as the observations, and analyzed these with the same 
analysis pipeline used for the actual observations. The simulated in-transit spectra are based 
on a "template" spectrum representing the disk-integrated spectrum of the star (described 
below), which we broaden to match the rotational velocity of HAT-P-lJfl We subtract from 
this template spectrum an unbroadened copy that is scaled by e and Doppler-shifted by v p , 
and then measure the radial velocity anomaly Av. We repeat this process for a grid of {e, v p } 
and approximate Av(e, v p ) with a two-dimensional polynomial fit. Differ ential rotation was 



ignored, as its effects are expected to be negligible (IGaudi fc Winnl 120071 ). 



3 The broadening kernel depends on the assumed limb-darkening law of the star, which we took to be a 
linear law with limb-darkening parameter u = 0.67. The results for the RM calibration formula are insensitive 
to the choice of u; very similar results were obtained for the choices u — 0.2 and u = 0.8. 
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The template spectrum should be similar to that of HAT-P-1 but with narrower lines 
because of the lack of rotat ional broadening ex hibited by the sub-planet spectrum. We 
selected the NSO solar atlas (IKurucz et al.lll984T) and a K eck/HIRES spectrum of HD 34411 
(T e// = 5911 K, [Fe/H] = +0.12: IValenti fc Fischer! hood ) . We found that the results based 



on either template are consiste nt with the functio n Av 



-ev p . This function is consistent 



with the analytic expressions of lOhta et al.l (120051 ) and iGimenezI (120061 ). We have found the 
best functional form of this "RM calibration" to vary from system to system; we also found 
a linear relation for TrES-2 (Winn et al. 2008). Thus, for this study, the calculated radial 
velocity of the star was taken to be the sum of the radial velocity of the Keplerian orbit, and 
the anomalous velocity Av = —ev p . 

The model had 12 free parameters: K, R*/a, Rp/R*, i, T c , v sini*, A, three additive 
constants for the 3 different groups of velocity data, and two limb- darkening coefficients az 
and u (to be explained in the next few paragraphs, see also Table H]). The orbital period 
was fixed at the value determined previously and the eccentricity was fixed at e = 0. The 
two model parameters relating to the RM effect are the line-of-sight stellar rotation velocity 
(vsinz*), and the angle between the projected stellar spin axis and orbit normal (A). The 
projected spin-orbit angle A is measured counterclockwise on the sky from the projected stel- 
lar rotational angular-momentum vector to the projected orbital angular-momentum vector 
(see Ohta et al. 2005 or Gaudi & Winn 2007 for a diagram). Due to the symmetry of the 
situation, a configuration with inclination i and spin-orbit angle A cannot be distinguished 
from a different configuration with inclination 180° — % and spin-orbit angle —A. To break 
this degeneracy we restrict i to the range from zero to 90 degrees, and allow A to range from 
-180° to +180°. 

For the photometric model, the limb-darkening law was assumed to be quadratic, I / J = 
l — a z (l — fi) — bz(l — ^) 2 , where /i is the cosine of the angle between the line of sight and the 
local surface normal. Given the precision of our data it is not possible to place meaningful 
constraints o n both az and bz- We fixed b z = 0.34, based on interpolation of the tables by 
Claretl (120041 ) . We allowed az to be a free parameter subject to a mild a priori constraint, 
shown in Eq. below, that enforces agreement with the tabulated value within «0.2. 
The choice of 0.2 is somewhat arbitrary and is fairly conservative, in the sense that better 
agreement is usually observed between fitted and theoretical limb darkening coefficients for 
the cases when such comparisons can be made (see, e.g., Winn, Holman, & Roussanova 2007 
and Southworth 2008). This approach is intermediate between the extreme approaches of 
fixing the limb darkening parameters exactly (placing too much trust in tabulated values) 
and allowing them to be completely free parameters (disregarding all theoretical knowledge 
of stellar atmospheres and possibly allowing unphysical parameter values). 
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For the RM model, we adopted a linear law [I/Iq = 1 — u(l — //)] for simplicity, since 
a quadratic law does not seem justified by the precision of the RM data. The appropriate 
choice of the limb-darkening coefficient u is not obvious. The Doppler-shift measurement is 
based on the portion of the spectr um between 50 00 and 6200 A, where the iodine absorption 
lines are plentiful. The tables of IClaretl (120041 ) lead to an expectation u ~ 0.67 for this 
spectral region. However, the Doppler information arises primarily from the steep sides of 
the stellar absorption lines in this region, and the degree of limb darkening in the lines may 
differ from the degree of limb darkening in the continuum, since the line radiation arises 
from a different depth in the stellar atmosphere. 

To investigate this issue we examined a Kurucz (1979) ATLAS12 plane-parallel model 
stellar atmosphere with T eS = 5750 K, \ogg = 4.5, [Fe/H] = 0.0)0, which was originally 
computed for the star XO-1 (McCullough et al. 2007) but whose properties are a reasonable 
match to those of HAT-P-1 (T cff = 5975 K, logg = 4.5, [Fe/H] = +0.1; Torres et al. 2008). 
The stellar intensity was computed for 17 different values of //, with a resolving power of 
500,000. When the spectrum is averaged betwe en 500-620 n m, the best-fitting linear limb- 
darkening coefficient is 0.65, in agreement with IClaretl (12004! ). According to the model, the 
degree of limb darkening is smaller in the steepest portions of absorption lines. For example, 
for a band centered on one of the Mg b triplet lines at 518.5 nm, we find u = 0.66 when 
the bandwidth is AA = 1 nm, and u = 0.50 when AA = 0.02 nm (encompassing only the 
steepest portion of the line). For other strong lines we also find that u is decreased by 0.1-0.2 
in the cores. For this reason, we chose to allow a to be a free parameter, with the same type 
of a priori constraint used for the photometric limb-darkening law (see below). 

The fitting statistic was 



X 



287 

E 

3=1 



fj(obs) - /j(calc)' 



125 



Vj(obs) — v j(calc) 



a 



v,3 



+ 



a z - 0.18 
(L2 



+ 



u 



0.67 



0.2 



where /^(obs) are the relative flux data from the composite light curve and afj is the out-of- 
transit rms. Likewise Vj (obs) and (7 y^j 21X6 the radial- velocity measurements and uncertainties 
after adding the jitter as described above. The last two terms represent a priori constraints 
on the linear limb-darkening coefficients az (for the photometric data) and u (for the radial- 
velocity data). As before, we solved for the model parameters and uncertainties using a 
Markov Chain Monte Carlo algorithm. We used a chain length of 10 6 steps and adjusted 
the perturbation size to yield an acceptance rate of ~40%. The posterior probability distri- 
butions for each parameter were roughly Gaussian, so we adopt the mean as the "best-fit" 



1) 



4 Downloadcd from kurucz . harvard . edu. 
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vable and the standard deviation as the 1-a error. For the joint model fit the minimum x 2 is 
412.2, with 402 degrees of freedom, giving xl = 1-03. This nearly "perfect" goodness-of-fit 
statistic should be interpreted as a check on the appropriateness of our data weights, rather 
than an independent check on the validity of the model, because we inflated the RV errors 
and attributed the RV scatter to jitter, and likewise we set the flux uncertainties equal to 
the out-of-transit RMS flux. 



4. Results 



The results from our analysis are given in Table |U The parameters depending on the 
RV data are fsinz* = 3.74 ± 0.30 km s _1 and A = 3?7 ± 2?1. (Below, we argue that the 
true error in vsini* is subject to an additional systematic error of 0.5 km s -1 .) The small A 
indicates close alignment between the sky-projected stellar spin axis and orbit normal. Our 
measured rotation velocity is higher than the value vsini* = 2.2 ± 0.2 km s _1 that was re- 
ported by lBakos et all (120071 ). Using the Spectroscopy Made Easy (SME) program (Valenti 
& Piskunov 1996, Valenti & Fischer 2005), we reanalyzed the same HIRES template obser- 
vation of HAT-P-1 that was analyzed by Bakos et al., and found v sin i± = 3.4 ± 0.5 km s _1 , 
which agrees with the result of our RM analysis. All other spectroscopic parameters from our 
SME analysis agreed with those reported previously. The primary difference in our analysis 
is that we used an appropriately narrow instrumental profile width, or equivalently a higher 
resolution, as measured from the model instrumental profile used in our Doppler analysis. 
The lower resolution assumed by Bakos et al. artificially compensated for rotational broad- 
ening, resulting in an erroneously low value of v sin (Debra Fischer private communication 
2008). 



The parameters t hat re ly pr imarily on our photo metry agree well with those of l Winn et al 



( )2007al ) , ISouthworthl (120081 ) and iTorres et all (120081 ) ; and those authors generally agree with 
one another, though there are differences in the exact treatment of red noise and limb 
darkening. As an additional check on our quoted errors, and in particular our assumption 
that the composite light curve had uncorrelated photometric errors, we applied the "residual- 
permutation" or "rosary-bead" method. In this method, one calculates the residuals between 
the photometric data and the best-fitting model, and creates many different "realizations" 
of the data that preserve any time-correlated noise by time-shifting the residuals and adding 
them back to the model. Then, the a posteriori probability distribution for each parameter 
is estimated by minimizing x 2 f° r each different realization of the data, and creating his- 
tograms of the parameter values. The error bars returned by this method were similar to, 
or smaller than, the error bars quoted in Table 4. 
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The result for the photometric limb-darkening parameter is az = 0.26 ± 0.06, showing 
that the data prefer a slightly more limb-darkened star than in the ATLAS models from 
which limb-darkening coefficients were tabulated by I Clare tl (120041 ) . The result for the radial- 
velocity limb-darkening parameter is u = 0.901q'2q. This is about la larger than the expected 
continuum value of 0.67, even though the ATLAS models predict that the steepest portion of 
the absorption lines (which provide most of the Doppler information) should exhibit smaller 
limb darkening. Assuming the models are correct, it is possible that the limb-darkening 
parameter u is compensating for an inaccuracy in our model of the RM effect, especially for 
the ingress and egress phases where limb darkening is strongest. This issue deserves further 
investigation, perhaps by increasing the sophistication of our RM calibration procedure (see 
§ I3.3p . using spatially resolved theoretical intensity distributions rather than an empirical 
stellar template. Fortunately this issue affects only the results for v sin i*, and not for A. This 
can be understood because u and v sin i+ both depend on the amplitude of the anomalous 
Doppler shift, while A depends almost entirely on the timing of the null of the anomalous 
Doppler shiftjj We verified this by fixing u at values between 0.5 and 0.9 and observing that 
the results for usinz* change by 0.5 km s _1 while the results for A are unchanged. Thus, 
we conclude that our result for usinz* is subject to a systematic error of approximately 
0.5 km s . In Table 4, we have added this systematic error in quadrature to the statistical 
error of 0.30 km s _1 giving a total error of 0.58 km s _1 . 



5. Summary and Discussion 

We have obtained high-precision photometric and spectroscopic measurements of the 
star HAT-P-1. Our in-transit spectroscopic observations clearly show the anomalous Doppler 
shift due to the Rossiter-McLaughlin effect, and we find that the angle between the sky 
projections of the stellar spin axis and the orbit normal is 3?7 ± 2?1. Additional Doppler 
measurements made during non-transit orbital phases allow us to constrain the orbital ec- 
centricity to e < 0.067 with 99% confidence. We measured the transit times from two new 
light curves and refined the orbital period by nearly an order of magnitude. 

The HAT-P-1 system is an interesting case for pla netary migration t heories because it 



is known to have a stellar companion in a wide orbit (IBakos et al.ll2007t ). a key ingredient 



for the Kozai mechanism. It is also suggestive that the radius of the planet HAT-P-lb is on 



5 For other transiting systems with small impact parameters, such as TrES-1 (Narita et al. 2007) and 
HAT-P-2 (Winn et al. 2007, Loeillet et al. 2008), there is a strong degeneracy between vsini* and A (Gaudi 
& Winn 2007). This type of systematic error may affect the results for A in those cases. 
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a relic of tidal energy dissipation (IFabrycky fc Tremaine 



the high end of theoretical expect ations (IBakos et al.l 120071; IWinn et al.ll2007al ). which may 



20071 ). Had the planetary system 



exhibited a large spin-orbit misalignment, it would have provided evidence for a scenario 
in which HAT-P-lb migrated to its current orbit a s a result of Kozai oscillatio ns, coupled 
with tidal dissipation within the planetary interior (IFabrycky fc Tremaind 120071 ). However, 
the small value of A does not necessarily rule out t he Kozai migration scenario, as small 
spin-orbit angles are not exc luded by the simulations (IFabrycky fc Tremaindl2007l ; IWu et al. 
20071 ; iNagasawa et al.l 120081 ) . But taken together with the 7 other planetary systems with 
small values of A, it seems likely that the dominant migration mechanism responsible for hot 
J up iters preserves spin-orbit alignment. 

Even if gravitational few-body mechanisms such as Kozai cycles do not represent the 
dominant migration channel for the formation of hot Jupiters, these mechanisms may be 
nonetheless be responsible for con figurations of s ome of close-in planets. A prime example is 
the orbit of HD 17156b, for which iNarita et al.l (l2007bl ) reported a 2.5a detection of a large 
misalignment (A = 62° ±25°). Since small values of A are not excluded by any of the existing 
theories of planet migration, misaligned systems like HD 17156 provide the most important 
tests of the various migration mechanisms. With wide-field transit surveys discovering new 
planets at an accelerating pace, the sample of measured spin-orbit angles will soon be large 
enough and precise enough to directly confront the theory. 
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Table 1. Doppler Shift Measurements of HAT-P-1 



Telescope Codc a Heliocentric Julian Date Radial Velocity Measurement Uncertainty 

(m s _1 ) (m s _1 ) 



1 


2453927.06954 


-50.33 


1.66 


1 


2453927.96670 


-41.12 


1.73 


1 


2453931.03823 


-23.27 


1.64 


1 


2453931.94201 


-57.75 


1.84 


1 


2453932.03731 


-55.85 


1.87 


1 


2453933.00132 


-9.75 


1.84 


1 


2453933.92609 


53.85 


2.32 


1 


2453934.90529 


38.18 


3.75 


1 


2453934.90853 


35.21 


2.09 


1 


2453960.99790 


62.44 


1.87 


1 


2453961.99776 


5.91 


1.83 


1 


2453963.09343 


-57.04 


2.09 


1 


2454287.93861 


9.85 


4.36 


1 


2454287.94541 


8.74 


2.75 


1 


2454287.95252 


13.18 


3.04 


1 


2454287.96032 


18.46 


2.68 


1 


2454287.96650 


2.63 


2.75 


1 


2454287.97330 


9.02 


2.81 


1 


2454287.97916 


-1.48 


4.36 


1 


2454287.98299 


8.69 


3.50 


1 


2454287.98721 


13.22 


3.54 


1 


2454287.99136 


20.22 


3.31 


1 


2454287.99485 


15.47 


3.78 


1 


2454287.99928 


15.27 


4.72 


1 


2454288.00334 


20.35 


4.62 


1 


2454288.00728 


22.05 


3.93 


1 


2454288.01142 


17.23 


3.70 


1 


2454288.01552 


15.48 


3.70 


1 


2454288.01886 


11.54 


3.33 


1 


2454288.02300 


23.50 


3.29 


1 


2454288.06101 


-13.13 


2.45 


1 


2454288.06483 


-33.06 


2.73 


1 


2454288.06840 


-27.90 


2.75 


1 


2454288.07171 


-19.28 


2.48 


1 


2454288.07475 


-11.57 


2.47 


1 


2454288.07773 


-11.96 


2.52 


1 


2454288.08068 


-14.12 


2.42 


1 


2454288.08364 


-21.86 


2.36 


1 


2454288.08669 


-6.65 


2.47 


1 


2454288.08978 


-1.64 


2.57 


1 


2454288.09297 


-0.91 


2.38 


1 


2454288.09591 


-4.42 


2.29 


1 


2454288.09896 


4.43 


2.54 


1 


2454288.10208 


-2.81 


2.57 


1 


2454288.10554 


1.99 


2.29 
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Table 1 — Continued 



Telescope Codc a Heliocentric Julian Date Radial Velocity Measurement Uncertainty 

(m s _1 ) (m s _1 ) 



1 


2454288.10893 


-0.73 


2.53 


1 


2454288.11884 


-5.13 


2.46 


1 


2454288.12233 


-1.83 


2.54 


1 


2454288.13167 


-10.34 


2.76 


1 


2454288.13551 


-10.33 


3.42 


1 


2454345.78243 


20.49 


2.11 


1 


2454345.78646 


20.04 


2.59 


1 


2454345.78999 


12.91 


2.58 


1 


2454345.83003 


13.35 


2.43 


1 


2454345.96340 


13.66 


2.66 


1 


2454346.02412 


8.30 


2.60 


1 


2454346.02708 


9.43 


2.52 


1 


2454346.02998 


6.72 


2.49 


1 


2454346.03951 


14.27 


2.59 


1 


2454346.04242 


18.06 


2.59 


1 


2454346.04529 


17.30 


2.89 


1 


2454346.04806 


18.58 


2.56 


1 


2454346.05089 


16.43 


2.69 


1 


2454346.05366 


12.71 


2.72 


1 


2454346.05642 


37.09 


2.57 


1 


2454346.05914 


16.46 


2.59 


1 


2454346.06189 


21.25 


2.52 


1 


2454346.07123 


5.91 


2.59 


1 


2454346.07399 


7.51 


2.61 


1 


2454346.07681 


15.06 


2.53 


1 


2454346.07969 


0.49 


2.43 


1 


2454346.08252 


0.13 


2.50 


1 


2454346.08557 


7.90 


2.51 


1 


2454346.08879 


-14.99 


2.54 


1 


2454346.09189 


-9.72 


2.68 


1 


2454346.09493 


-3.40 


2.79 


1 


2454346.09819 


-9.89 


2.73 


1 


2454346.10137 


-17.76 


2.61 


1 


2454346.10442 


-18.61 


2.72 


1 


2454346.10748 


-12.82 


2.76 


1 


2454346.11065 


-20.48 


2.55 


1 


2454346.11385 


-24.27 


2.52 


1 


2454346.11714 


-24.44 


2.62 


1 


2454346.12030 


-25.08 


2.61 


1 


2454346.12342 


-21.52 


2.45 


1 


2454346.12652 


-22.09 


2.61 


1 


2454346.12959 


-23.63 


2.63 


1 


2454346.14545 


-2.77 


2.79 


2 


2454318.07432 


2.81 


7.04 


2 


2454318.08202 


17.62 


7.33 
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Table 1 — Continued 



Telescope Code a 


Heliocentric Julian Date 


Radial Velocity 

(ms- 1 ) 


Measurement Uncertainty 
(ms- 1 ) 


2 


2454318.11099 


24.15 


6.57 


2 


2454349.00272 


6.64 


5.71 


2 


2454349.01034 


13.10 


6.62 


2 


2454349.01795 


2.20 


4.55 


2 


2454349.02557 


5.69 


5.95 


2 


2454349.03319 


9.39 


4.64 


2 


2454363.77654 


-40.68 


7.64 


2 


2454363.78414 


-37.33 


6.61 


2 


2454363.90169 


-32.80 


6.31 


2 


2454363.90931 


-25.78 


7.03 


2 


2454363.91692 


-25.12 


6.14 


2 


2454363.92455 


-33.95 


8.05 


2 


2454363.93216 


-37.52 


7.64 


2 


2454363.93978 


-49.54 


7.06 


2 


2454363.94738 


-48.82 


6.71 


2 


2454363.95500 


-58.17 


7.52 


2 


2454363.96261 


-49.23 


6.98 


2 


2454363.97023 


—66.33 


7.15 


2 


2454363.97785 


-65.61 


5.92 


2 


2454363.98547 


-64.40 


7.38 


2 


2454363.99307 


-61.59 


7.59 


2 


2454364.00069 


-62.49 


6.41 


2 


2454364.00830 


-50.84 


7.51 


2 


2454364.01591 


-44.29 


7.59 


2 


2454364.03113 


-47.85 


9.03 


2 


2454364.03874 


-53.34 


8.92 


2 


2454364.04636 


-56.22 


7.92 


2 


2454364.05397 


-62.86 


7.99 


2 


2454364.06158 


-55.24 


8.02 


2 


2454364.07680 


-56.72 


8.92 


2 


2454364.08441 


-52.38 


8.53 


3 


2453897.11172 


-2.89 


5.00 


3 


2453899.12580 


52.53 


4.77 


3 


2453900.11917 


-34.51 


4.65 


3 


2453901.10176 


-37.51 


4.88 



a (l) HIRES, Keck I 10m telescope, Mauna Kea, Hawaii. (2) HDS, Subaru 8m telescope, Mauna 
Kea, Hawaii, based on observations and data reduction procedures described in this work. (3) HDS, 
Subaru 8m telescope, Mauna Kea, Hawaii, based on observations and data reduction procedures 
described by Bakos et al. (2007). 
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Table 2. Photometric Measurements of HAT-P-1 



Telescope Code b Heliocentric Julian Date Relative Intensity 



1 


2454381.71537 


0.9988 


1 


2454381.71588 


1.0002 


1 


2454381.71640 


1.0007 


1 


2454381.71691 


1.0028 


1 


2454381.71904 


1.0018 


1 


2454381.71955 


1.0002 



a Full table available in online addition. 

b (l) Nickel telescope, Mt. Hamilton, California. (2) MAGNUM 
2m telescope, Haleakala, Hawaii. 
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Table 3. Midtransit times of HAT-P-1 



Midtransit time [HJD] Reference 



2453984. 


.39700 ± 


0. 


00900 


1 


2453979. 


.92994 ± 


0. 


.00069 


2 a 


2453988. 


.86197 ± 


0. 


.00076 


2 


2453997. 


.79200 ± 


0. 


.00054 


2 


2453997. 


.79348 ± 


0. 


.00047 


2 


2454006. 


72326 ± 


0. 


00059 


2 


2454015. 


.65338 ± 


0. 


.00107 


2 a 


2454069, 


.23795 ± 


0. 


.00290 


2 


2454363. 


.94601 ± 


0. 


.00091 


3 


2454381. 


.80849 ± 


0. 


.00125 


3 



Note. - References: (1) Bakos et 
al. (2007), (2) Winn et al. (2007), (3) 
This work. 

a The midtransit time for this event 
that was originally reported by Winn et 
al. (2007) was incorrect. The data re- 
ported here have been corrected. 
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Table 4. System Parameters of HAT-P-1 



Parameter 



Value 



Uncertainty 



Transit ephemeris: 
Orbital period, P [d] 
Midtransit time [HJD] 



4.4652934 
2454363.94656 



0.0000093 
0.00072 



Rossiter- McLaughlin parameters: 

Projected spin-orbit angle, A [deg] 

Projected stellar rotation rate, vsini* [km s -1 ] 



3.7 
3.75 



2.1 
0.58 



Photometric transit parameters: 

Planet-to-star radius ratio, R p / R* 

Orbital inclination, i [deg] 

Scaled semimajor axis, a/R* 

Transit impact parameter 

Transit duration [hr] 

Transit ingress or egress duration [hr] 



0.11295 
86.28 
10.67 
0.693 
2.798 
0.510 



0.00073 
0.20 
0.25 
0.023 
0.019 
0.031 



Spectroscopic orbital parameters: 

Orbital eccentricity, e 

Velocity semiamplitude, K [m s _1 ] 

Additive velocity (Keck/HIRES) [m s" 1 ] 

Additive velocity (Subaru, this work) [m s _1 ] 

Additive velocity (Subaru, Bakos et al. 2007) [m s _1 ] 




59.3 
1.26 
-45.3 
13.0 



assumed b 
1.4 
0.66 
1.5 
3.1 



Derived system parameters: 
[M ] b 

R* [R ] c 
M p [M Jup ] c 

R P [Pju P ]° 

M p /M, c 



1.133 
1.115 
0.524 
1.225 
0.000441 



0.077 
0.050 
0.031 
0.059 
0.000020 
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a When the orbital eccentricity is taken to be a free parameter, we find an upper limit 
of 0.067 with 99% confidence. 

b The stellar mass was not determined from our data. The value given here is from 
Torres, Winn, & Holman (2008). 

c This result depends on the assumed value and uncertainty in the stellar mass. 
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Fig. 1. — Top: Relative radial velocity measurements of HAT-P-1, from this work and 
from Bakos et al. (2007). The solid line is the best-fitting model. The typical measurement 
uncertainties are illustrated as points with error bars in the lower left corner. A detailed view 
near midtransit is shown in Figure HI Bottom: Residual radial velocities after subtracting 
the best-fitting model. 
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Fig. 2. — Photometry of transits of HAT-P-1, using the Nickel lm telescope and a Z-band 
filter (top), and the MAGNUM 2m telescope and a V-band filter (bottom). The solid lines 
show the best-fitting model. 
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Fig. 3. — Residuals from a linear ephemeris that was fitted to the transit times in Table [3j 
The best-fitting ephemeris has a period P = 4.4652934 ± 0.0000093 days and midtransit 
time T c = 2454363.94656 ± 0.00072 (HJD). 
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Fig. 4. — Top: Relative radial velocity measurements of HAT-P-1, centered on the midtran- 
sit time. Bottom: Composite, time-binned transit light curve of HAT-P-1, based on the Z 
and z-band data from Winn et al. (2007) and this work. The rms scatter of the residuals is 
0.00057. This composite light curve was fitted simultaneously with the radial-velocity data. 



